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Abstract 

Sensory hair cells in auditory and vestibular organs rely on active mechanisms to achieve high sensitivity 
and frequency selectivity. Recent experimental studies have documented self-sustained oscillations in hair 
cells of lower vertebrates on two distinct levels. First, the hair bundle can undergo spontaneous mechani- 
cal oscillations. Second, somatic electric voltage oscillations across the baso-lateral membrane of the hair 
cell have been observed. We develop a biophysical model of the bullfrog's saccular hair cell consisting of 
two compartments, mechanical and electrical, to study how the mechanical and the voltage oscillations 
interact to produce coherent self-sustained oscillations and how this interaction contributes to the overall 
sensitivity and selectivity of the hair cell. The model incorporates nonlinear mechanical stochastic hair 
bundle system coupled bi-directionally to a Hodgkin-Huxley type system describing somatic ionic cur- 
rents. We isolate regions of coherent spontaneous oscillations in the parameter space of the model and 
then study how coupling between compartments affect sensitivity of the hair cell to external mechanical 
perturbations. We show that spontaneous electrical oscillations may enhance significantly the sensitivity 
and selectivity of the hair cell. 



Introduction 

Hair cells are peripheral sensory receptors in the inner ear of vertebrates which transduce mechanical 
vibrations to electrical signals [1J2] . External mechanical perturbations are translated to electrical signals 
through the mechano-electrical transduction (MET) channels located in the stereocilia, apical projections 
of the cell immersed in potassium (K) rich endolymph. Stereocilia are arranged in the hair bundle in 
rows of increasing heights and linked by the so-called tip links. External mechanical perturbations deflect 
the hair bundle resulting in changing of the tension in the tip links, which may elicit opening or closing 
of MET channels, mediating the influx of potassium and calcium ions to the cell and thus changing 
the cell's potential. Experimentally observed extreme sensitivity and frequency selectivity of hair cells 
suggested several mechanisms of active amplifications [l}j3]. One of these mechanisms relies on active 
dynamics of the hair bundle which can be placed at the verge of dynamical instability, i.e. Andronov- 
Hopf bifurcation, providing the cell with high sensitivity, sharp frequency tuning curve and a compressive 
nonlinearity [4j|6]. MET channels in stereocilia opening in concert exert forces on the hair bundle, the 
phenomenon known as gating compliance [7], resulting in a nonlinear force-displacement relation. In 
some hair cells the gating compliance may result in a differential negative stiffness of the hair bundle 
leading to its instability |8 . Combined with myosin-mediated adaptation [9 10 this nonlinearity result 



in spontaneous mechanical hair bundle oscillatory motion. In particular, spontaneous oscillations of the 
hair bundle have been observed in many lower vertebrates such as amphibian sacculus, or turtle basilar 
paplia 11} 15 and presumably responsible for the phenomenon of otoacoustic emissions in nonmammals 
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161. Importantly, spontaneous hair bundle oscillations are essentially stochastic due to several 



thermal noise sources 17 . This natural stochasticity can be reduced in coupled hair bundles 18 19 . 

Intriguingly, the same type of hair cells show spontaneous oscillations of the membrane potential. 
For example, several experimental groups reported on spontaneous voltage oscillations in saccular hair 
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cells 20-23 . The functional role of these large-amplitude voltage oscillations is currently not known. One 
suggestion is that voltage oscillations provide additional feedback to the hair bundle dynamics which may 
help to reduce thermal fluctuations and thus to improve sensitivity of the hair cell 



24 . Indeed earlier 



studies reported on reverse electro-mechanical transduction, providing evidence for bi-directional coupling 
between the hair bundle complex and the potential of the cell bo dy [25] . Furthermore, transepithelial 
electrical stimulation of hair cells elicits the hair bundle motion [26] . Finally, a recent experimental 
study showed that somatic ionic conductances affect significantly the mechanical oscillations of hair 
bundles. 
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In this paper we develop a model of sensory hair cell which includes mechanical and electrical compart- 
ments coupled bidirectionally. Sensitivity, frequency selectivity and compressive nonlinearity are three 
main features of the hair bundle dynamics [6] . We use the model to study how the coupling affects these 
three feature of the hair cell dynamics. 



Two compart mental model of a hair cell 

Our two compartmental system is based on previously developed models for the hair bundle mechanics 
and for the dynamics of the cell potential of the bullfrog's saccular hair cells. The coupling between the 
two compartments is bidirectional. The coupling from the mechanical to electrical compartment is due 
to direct mechano-electrical transduction, i.e. due to the MET current [28|[29| . The coupling from the 
electrical to mechanical compartment is mediated by the reverse electro-mechanical transduction whereby 
the cell voltage affects the calcium ions concentration near the myosin molecular motors controlling the 
rate of adaptation of the hair bundle. 

Several simple models were proposed to account for the active hair bundle motion including spon- 
taneous mechanical oscillations in the bullforg's sacculus 14 17[ 30 33 . In this paper we use a model 



proposed in |T7[[3TJ . The overdamped dynamics of the hair bundle is described by two stochastic differ- 
ential equations: first for the position of the bundle tip, X and second for the displacement of the myosin 
motors, X a along the stereocilia: 

= -K GS (X -X a -DP Q )- K SP X + F ext (t) + ^2k B T\ v (t), (1) 

dt 

A a ^ = K GS (X -X a -DP )- F max (l -SP ) + v/3fc B TA a ?7 a (t), 
dt 

p (Y X ) - 1 A - P l^G+K GS D 2 /(2N)]/(k B T) 

o[ ' a> i + Ae-( x - x *) K GsD/(Nk B T)> 

In Eqs.Q rj{t) and T] a (i) are uncorrelated Gaussian white noise sources modeling thermal noise due to 
Brownian motion, MET channel cluttering and stochastic binding and unbinding of adaptation motors 
|17| . The external stimulus enters the model via the external force term F cxt (t) . Definition of parameters 
and their values are given in Table [I] and are essentially the same as used in 17 18 . The control 



parameters of the hair bundle model are 5" and F mra which determine the strength of calcium (Ca) 
influence on adaptation motors and maximal force generated by the motors at stall when all MET 
channels are closed, respectively. 

The dynamics of the cell's potential can be described by a Hodgkin-Huxley type system developed 



21-23 . We used a modified version of this model described in details in [34 . The model includes 6 



basolateral ionic currents and the MET current. The current balance equation is 

C m — = — JkI — ^BKS — ^BKT — ^DRK — — ^Ca — ^MET, (2) 

dt 

where Jki is the inwardly rectifier potassium (K) current, Ibks & n d Ibkt are Ca-activated steady and 
transient K currents, 1^ is sodium / potassium h-type current, Idrk is the direct rectifier K current and 
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Table 1: Parameters of the mechanical compartment 



Symbol 


Definition 


value 


A 


Drag coefficient of the hair bundle 


2.8 /xNs/m 


A a 


Effective friction coefficient of the molecular motors 


fO.O /xNs/m 


Kgs 


Gating spring stiffness 


0.75 mN/m 


K SP 


Stereocilia pivots stiffness 


0.6 mN/m 


D 


Gating swing of MET channel 


60.9 nm 


T 


Ambient Temperature 


300K 


± max 


Maximal motors force 


55 pN 


s 


Dimensionless Ca feedback strength 


1.13 


Po 


Open probability of MET channels 




N 


Number of MET channels 


50 


AG 


Energy change on MET channel opening 


10k B T 



lea, is the Ca current. The equation for the potential ([2| is accompanied by equations for kinetics of ionic 
currents listed and for intracellular [Ca 2+ ], totaling 12 differential equations. The detailed description 
of this system and parameters is provided in [34] . The control parameters of the electrical compartment 
are the maximal conductance of inwardly rectifier current, 3ki, and relative strength of Ca-activated 
currents, &k- Depending on these two parameters the electrical compartment shows rich variety of 



oscillatory patterns including quasiperiodic, bursting and chaotic oscillations 34 . By setting the fast 
activation variable for the inwardly rectifier current (Iki) and [Ca 2+ ] to their steady states we reduced 
the dimension of the electrical compartment by 2. We verified that this reduction did not change the 
voltage dynamics quantitatively. 

Forward coupling from the mechanical to the electrical compartment is due to the MET current 
entering in Eq.|2]), 

^MET = gMETPo(X, X a )(V(t) — -EmEt), (3) 

where £"met = mV is the reversal potential potential and <?met is the maximal value of the MET 
conductance, i.e. when all MET channels are open. The MET conductance ^met thus serves as a 
"forward" coupling parameter. The dynamics of the hair bundle depends crucially on the concentration 
of Ca ions at the adaptation motors cite [Ca 2+ ]M, which is determined by the open probability of MET 



channels and by the potential of the cell V via electrodiffusion 14 31 . For example, a decrease of 
the receptor potential V enhances the influx of Ca ions through MET channels resulting in increase of 
[Ca 2+ ]M which in turn inhibits the motor activity. This backward 'voltage to mechanics" coupling can 
be introduced via voltage-dependent calcium feedback parameter S(V) in Eqs.Q [24] , 

S W = S °J(yly f{V 1 = T^V> e = 2q e /(k B T), 

where Sq is a value of the feedback parameter S at a reference potential V — Vq — — 55 mV and q e is the 
elementary charge. Since the voltage deviations are in the range -80 - -30 mV, S(V) can is approximated 
by its first order Taylor expansion around the reference potential, 

S(V) =S a [l + a(V/V - 1)] , (4) 

where we introduced the dimensionless parameter a controlling the strength of the backward coupling, 
such that a = corresponds to the hair bundle uncoupled from the variations of the receptor potential. 

Two compartmental stochastic model composed of the hair bundle model ([lj and of the Hodgkin- 
Huxley type model for the membrane potential ^ coupled via Eqs.(|3]) and Q was studied numerically 
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using the Eulcr-Maruyama method with the time step 0.1 ms. The control parameters ol both com- 
partments were chosen such that in the absence of coupling (a = 0, <7met — 0), stimulus and noise, the 
hair bundle and the voltage were at a stable equilibrium. In particular we fixed the parameters of the 
mechanical compartment at (F max , S) = (55 pN, 1.13), and the parameters of the electrical compartment 
at (&Ki <?Kl) — (0.01, 1.0 nS). The coupling strengths (?met and a were then used as control parameters. 
Bifurcation analysis of equilibrium states of the deterministic model, i.e. stochastic terms rj, r) a = in 
Eqs{l] was performed using CONTENT continuation software package 35 . We then explored sponta- 



neous stochastic dynamics of the model versus the coupling strengths. In particular we calculated the 
power spectrald density (PSD) of spontaneous hair bundle motion and then measured the quality factor 
Q of the main peak in the specrum and the mean amplitude of spontaneous oscillations as a function of 
the coupling strengths (?met and a. 

The sensitivity and frequency selectivity of the hair cell was estimated using periodic external force, 

F ext (t) = F cos(27r/ s t), (5) 

applied to the mechanical compartment ([I]). We calculated the time dependent averages of the hair bundle 
position, (X(t)) and of the membrane potential (V(t)} by averaging over 300 realizations of stochastic 
terms rj and rj a during 500 periods of the external sinusoidal force. We then estimated the mechanican 
and electrical sensitivities as the ratios of the first Fourier harmonic of these averages, X(f s ); V(f s ), to 
the amplitude of the stimulus, F 0l 

Xf(f s ,F Q ) = \X(f s )\/F Q , Xv(f s ,F ) = \V(f s )\/F , (6) 

in units of nm/pN for the mechanical sensitivity Xfi an d mV/pN for the electrical sensitivity, xv ■ 



Results and discussion 

In absence of coupling and noise the deterministic model is at a stable equilibrium, i.e. the mechanical 
compartment is at equilibrium with most of MET channels closed, P f=s 0.19 and the electrical compart- 
ment is at V = —53.5 mV. With the increase of coupling the equilibrium bifurcates to a limit cycle via 
the Andronov-Hopf (AH) bifurcation. On the parameter plane (^met, ot) the lines of the AH bifurcations 
isolate a region of spontaneous oscillations, shown in Fig. [T] In the region below the lower AH bifurcation 
line, most of MET channels are closed, while in the region above the upper AH bifurcation line most 
of MET channels are open and the cell is depolarized. In the region between the AH bifurcation lines 
the deterministic model exhibit large-amplitude synchronous oscillations in mechanical and electrical 
compartments. 

Thermal noise induces stochastic hair bundle oscillations even outside the deterministic oscillatory region. 
Figure [T] (left panel) shows that a region of large amplitude mechanical oscillations extends below the 
lower boundary of deterministic oscillatory region. Noise-induced bundle motion leads to opening of MET 
channels and since the Ca feedback is relatively weak below the lower AH bifurcation line, adaptation 
brings the bundle back to the equilibrium position allowing large-amplitude oscillation. On the contrary, 
above the upper boundary of deterministic oscillatory region large values of the backward coupling a lead 
to weak adaptation resulting in small-amplitude fluctuations of the hair bundle around the equilibrium 
with open MET channels. This is illustrated in Fig. [2] which shows the hair bundle displacement and 
the voltage traces at three points corresponding to three distinct regions on the parameter plane of the 
system (cf Fig. [I]). 

The coherence of spontaneous oscillations quantified with the quality factor, Q, of the main peak in the 
PSD of the hair bundle position, shows non-monotonous behavior being maximal at the center of the 
oscillatory region (point b in the right panel of Fig. nl). In the region below the lower AH line (point a 
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Figure 1: Spontaneous dynamics of the hair cell model versus coupling strengths between mechanical and 
electrical compartments. Left: Standard deviation (SD) of the hair bundle position, the color bar shows 
SD in nm. Right: Quality factor Q of the main peak in the power spectral density of the hair bundle 
oscillation. On both panels the AH bifurcation lines of the deterministic system are shown by white lines. 
Points a, b and c indicate three different regions discussed in the text. 
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Figure 2: Time series (left ) and PSDs (right) of spontaneous mechanical oscillations(red,[nm]), and 
membrane potential (black, [mV]) for the points shown in Fig. [l] Units of PSDs are [nm 2 /Hz] for the 
hair bundle displacement and [mV 2 /Hz] for the membrane potential. 



in Figs. [JJ and [2]) the dynamics of the system is dominated by noise resulting in noisy oscillations around 
the equilibrium with broad peaks in PSDs at the natural frequency determined by the imaginary parts 
of the equilibrium's eigenvalues. In the deterministic oscillatory region (point b in Figs, [JJ and [2]) both 
compartments show synchronous and coherent oscillations with sharp peaks in PSDs. Finally, the region 
above the upper AH line (point c in Figs. IT] and pi) is characterized by overdamped small amplitude 
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fluctuations with a low- frequency Lorentzian type PSD. The effect of the backward electro-mechanical 
coupling on the hair bundle dynamics is further illustrated in Fig. [3] showing that the coherence of 
spontaneous hair bundle oscillations is maximized by the backward coupling a. 
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Figure 3: Influence of the electrical compartment on the hair bundle oscillation. Left: Quality factor 
of the mechanical oscillation versus the backward coupling strength a for the fixed value of the MET 
conductance, <?met = 1 n S. Right: PSD of the hair bundle displacement for indicated values of a. 

Frequency response of mechanical and electrical compartments is shown in Fig. [4] for a weak mechanical 
sinusoidal stimulus Sensitivity of mechanical and electrical compartments peaks at the frequency of 
spontaneous oscillations. Furthermore, the peak values of sensitivity is largest for the coupling strength 
corresponding to the most coherent spontaneous oscillations (cf Fig. [31 left panel). The sensitivity as 
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Figure 4: Effect of the backward coupling on sensitivity of the hair cell model to external periodic force. 
Left column shows the mechanical sensitivity, i.e. the hair bundle sensitivity, \f- Right column shows the 
electrical sensitivity, \v- Top panels shows sensitivity versus stimulus frequency, / s , for the fixed value of 
the stimulus amplitude Fq = 2 pN and for indicated values of the backward coupling strength a. Bottom 
panels show the sensitivity versus stimulus amplitude, F . For indicated value of a the stimulus frequency 
was chosen to match the frequency of spontaneous oscillation. Direct mechano-electrical coupling strength 
was fixed at <?met = 1 nS. 

a function of the stimulus strength shown in the lower panels of Fig. [4j demonstrates the phenomenon 
of compressive nonlinearity. Both compartments respond linearly to weak stimuli, while the sensitivity 
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Figure 5: Maximal sensitivity of the hair cell model to weak sinusoidal force, Fq = 2 pN. Left panel shows 
maximal sensitivity of the hair bundle versus both forward (<7met) and backward (a) coupling strengths. 
Right panel shows sensitivities of the hair bundle and membrane potential versus a for the fixed value 
ffMET = 1 nS. 



to stronger stimuli is suppressed. The differences in response properties appear for strong stimuli Fq > 
100 pN where the hair bundle start to respond linearly again, while the sensitivity of the electrical 
compartment continues to decline. We note that similar effect was reported in 34 for voltage response 
where a linear model for the hair bundle was used. We also note that the backward coupling strength 
does not affect significantly the scaling of sensitivity with the stimulus amplitude. 

The hair bundle sensitivity was studied for a wide range of forward and backward coupling strengths 
(ffMET € [0, 3]nS), (a G [0, 2]). For each values of (/met and a we estimated the maximum sensitivity by 
varying the stimulus frequency / s and then plotted this maximum value of sensitivity as a color-coded 
map. The resulting maximal sensitivity map is shown in Fig. [5] left panel. This map clearly shows the 
existence of optimal coupling strengths at which the hair cell is most sensitive to periodic perturbations. 
For a fixed value of the forward coupling strength ,9met = 1-0 nS the sensitivity of the hair bundle peaks 
near a = 1.5, shown on the right panel of Fig.[5j and corresponding well to the most coherent spontaneous 
oscillation (see Fig. [3]) . 



Conclusions 

The role of active processes in operational performance of sensory hair cells is a topic of intense current 
interest in sensory neuroscience [T}|3| [36{[37] . In particular, the role of active hair bundle dynamics is a 
matter of debate [3j[6]- However, hair cells in lower vertebrates demonstrate also spontaneous oscillations 
of somatic potential which presumably may affect mechanics of the hair bundle, or may be a result of 
bi-directional coupling of mechanical and electrical constituents of the hair cell. Unlike many previous 
modeling works which studied mechanical and electrical dynamics separately, we constructed a model 
incorporating nonlinear mechanical and electrical compartments coupled bi-directionally. We showed that 
spontaneous oscillations may arise due to bi-directional coupling even when uncoupled compartments 
are quiescent. The coherence of oscillations can be enhanced by tuning the coupling strengths between 
compartments resulting in enhanced sensitivity and sharper tuning to weak periodic mechanical stimulus. 
At the same time, the model shows other distinctive behaviors, such as compressive nonlinearity. 
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